Accepted to ApJ (Version January 25, 2013) 

Preprint typeset using I^T^X style cmulatcapj v. 08/22/09 



X-RAY POINT SOURCE POPULATIONS CONSTITUTING THE GALACTIC RIDGE X-RAY EMISSION 

KUMIKO MORIHANA 1 , MASAHIRO TSUJIMOTO 2 , TESSEI YOSHIDA 3 , AND KEN EbISAWA 2,4 
Accepted to ApJ (Version January 25, 2013) 



o 

(N 
C 
>—> 

(N 



X 

43 

Oh: 

6 



> 

in 
oo 

O 

m 



ABSTRACT 

Apparently diffuse X-ray emission has been known to exist along the central quarter of the Galactic 
Plane since the beginning of the X-ray astronomy, which is referred to as the Galactic Ridge X-ray 
emission (GRXE). Recent deep X-ray observations have shown that numerous X-ray point sources 
account for a large fraction of the GRXE in the hard band (2-8 keV). However, the nature of these 
sources is poorly unders t ood. Using the deepest X-ray observations made in the Chandra Bulge Field 
(jRevnivtsev et all [2009, 20111) . we present the result of a coherent photometric and spectroscopic 
analysis of individual X-ray point sources for the purpose of constraining their nature and deriving 
their fractional contributions to the hard band continuum and Fe K line emission of the GRXE. Based 
on the X-ray color-color diagram, we divided the point sources into three groups: A (hard), B (soft and 
broad spectrum), and C (soft and peaked spectrum). The group A sources are further decomposed 
spectrally into thermal and non-thermal sources with different fractions in different flux ranges. From 
their X-ray properties, we speculate that the group A non-thermal sources are mostly AGNs and the 
thermal sources are mostly white dwarf (WD) binaries such as magnetic and non-magnetic cataclysmic 
variables (CVs), pre-CVs, and symbiotic stars, whereas the group B and C sources are X-ray active 
stars in flares and quiescence, respectively. In the logiV-logS curve of the 2-8 keV band, the group 
A non-thermal sources is dominant above wlO -14 ergs cm -2 s _1 , which is gradually taken over by 
Galactic sources in the fainter flux ranges. The Fe Ka emission is mostly from the group A thermal 
(WD binaries) and the group B (X-ray active stars) sources. 
Subject headings: Galaxy: bulge — Galaxy: disk — X-rays: diffuse background 



1. INTRODUCTION 

Since the dawn of the X-ray astronomy, the apparently 
diffuse emission of low surface brightness has been known 
to exist along the Galactic Plane (GP; \l\ <45°, \b\ <1°), 
which is referred to as the Galactic Ri d ge X-ray emis- 
sion (GRXE; e.g.. jWorrall et all [1981 [Warwick et all 
1985> IKovama et al.lll986ft . The X-ray spectrum is char- 
acterized by hard continuum with a strong Fe K emis- 
sion feature at 6-7 keV band (|Kovama et aLl figSS) . The 
GRXE is considered to have different origin s for the soft 
(<2 keV) and hard (>2 keV) emission (e.g.. lHands et al.l 
2004), and we focus on the la t ter in this paper. 

Recently, iRevnivtsev et al.l ([20091 ) has shown that the 
majority (~80%) of the Fe K-band emission was re- 
solved into point sources using the deepest X-ray ob- 
servations made with the Chandra X-ray Observatory 
( Wcisskopf ct al. 2002) at a slightly off-plane region of (I, 
6) = (0?113, -1?424) in the Galactic bulge {Chandra bulge 
field; CBF). The sensitivity of the study reached -10" 16 
ergs cm -2 s _1 in the hard (2-8 keV) X-ray band with a 
total integration time of nearly 1 Ms, which is >10 times 
deeper than all previous studies. In previous shallower 
surveys, only a small fraction of the GRXE was resolved 
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into point sources. For example, lEbisawa et al.l (2001, 
2005) conducted 100 ks Chandra observations at two 
partially overlapping regions in the GP at (I, 6)=(28?5, 
0?0), and found that only -10% of the total GRXE 
emission was resolved into point sources. These find- 
ings suggest that a large contribution is made for the 
GRXE by point sources in the flux range around 10~ 15 - 
10~ 16 ergs cm~ 2 s~ x . 

The nature of these numerous dim point sources re- 
mains unknown. They may be explained by the extrap- 
olation of the X-ray source population known in brighter 
flux ranges, or new classes of X-ray sourc es are required 
as m ajor contributers at this flux range (Ebisawa et al. 
2005). Understanding the X-ray source population at 
a certain flux range boils down to constructing the 
log N— log S curve separately for major classes of sources. 
In the hard X-ray band, the log N — log S curve was stud- 
ied with th e Advanced Satellite for Cosmology and As- 
trophysics (Sugizaki ct al. 2001) at a flux range brighter 
than ~10~ 125 ergs cm~ 2 s _1 an d with the Chandra 
and XMM-Newton Observatories (lEbisawa et al.l 120051 : 
lHands et aT]|2004t iMotch et al.ll2010D at a flux range be- 
tween ~10 -12 5 and ^10~ 145 ergs cm -2 s _1 . These stud- 
ies divided detected point sources phenomenologically 
based on the X-ray properties aided by optical and near- 
infrared identifications, and revealed that several classes 
of sources contribute to the log N — log S curve with dif- 
ferent fractions in different flux ranges. The brightest 
end of the curve is dominated by low-mass X-ray bina- 
ries, which saturates at ~10 -13 ergs cm -2 s _1 for being 
integrated all throughout our Galaxy in the line of sight. 
Cataclysmic variables (CVs) and X-ray active binary 
stars emerge as dominant contributers below the flux. 
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Extra-galactic sources, which are mostly active galactic 
nuclei (AGNs), co-exist with these Galactic sources in 
survey fields in all flux ranges. 

Even fainter flux range is accessible o nly with the deep- 
est Ch andra observation in the CBF. iRevnivtsev et all 
(|2011[ ) constructed the log N - log S curve with all point 
sources in the central 2f56 circle within 17f4xl7f4 fields 
(only 7% in the area) down to 10 -16 ergs cm~ 2 s _1 and 
discussed that the curve can be reproduced by the sum of 
CVs and active binary stars using the expected number 
of these sources in the line of sight and their luminos- 
ity fu nction derived in the Solar vicinity (jSazonov et al.l 
I2006D. However, in the prev ious studies in the CBF 
(|Revnivtsev et all 120091 120111 ), individual sources were 
treated collectively and their different X-ray properties 
were not considered, which should be useful to classify 
these sources. 

What we aim to do in this paper is to derive the X- 
ray photometric properties (flux and color) from all the 
sources in the entire CBF field and the spectroscopic and 
timing properties from bright sources through coherent 
data analysis. By using the entire field, we can increase 
the number of sources at a compensation of slightly lower 
averaged sensitivity in comparis on to the previous stud- 
ies on ly using the central part (Re vnivtsev et al.l 120091 
120111) . Then, we divide the sources phenomenologically 
into several groups based on these X-ray properties to 
construct the log TV — log S curve separately for each 
group to find their contributions at the flux range of in- 
terest. We also consider the contribution of these groups 
to the GRXE separately for the Fe K emission and the 
continuum emission in the hard band, as different classes 
are expected to have different equivalent width (EW) of 
the feature. 

The outline of this paper is as follows. In § [21 we 
present observations and data reduction. In § [3j we ex- 
tract point sources, derive the survey completeness limit, 
and obtain X-ray properties of individual sources. In §01 
we divide all the sources based on their X-ray colors and 
flux, and discuss their likely classes based on the spectro- 
scopic and timing properties of the brightest members of 
each group, and the composite spectrum of all sources in 
the group. We derive the contribution of these sources 
to the GRXE. We summarize the results in §0 For clar- 
ity, we use "groups" for subsets of detected point sources 
divided phenomenologically based on their X-ray prop- 
erties, and "classes" for different types of astronomical 
objects (e.g., AGNs, CVs, stars, etc) in this paper. 

2. OBSERVATIONS AND DATA REDUCTION 

We retrieved ten archived data of the CBF taken 
with the Advanced CCD Imaging Spectrometer (ACIS; 
iGarmire et alll2003h -l array on board Chandra. The ob- 
servations were carried out from 2008 May to August 
with a total exposure time of ~900 ks (Table [lj . The 
ACIS-I covers the 0.5-8.0 keV energy band with a spec- 
tral resolution of ^280 eV for the full width at a half 
maximum at 5.9 keV as of the observation dates 5 . The 
CCDs were operated with a frame time of 3.2 s. The data 
were down-linked with the very faint telemetry format. 

We retrieved pipeline products and reduced the data 

5 See http://cxc.harvard.edu/cal/Acis/ for details. 



TABLE 1 
Observation Log 



ObsID 


Date 


Coordinate (J2000.0) 


t a 

<<cxp 


Roll 






R.A. 


Decl. 


(KS) 




9500 


2008-07-20 


17:54:38.5 


-29:35:47 


162.56 


279.99 


9501 


2008-07-23 


17:54:39.1 


-29:35:53 


131.01 


278.91 


9502 


2008-07-17 


17:54:39.8 


-29:35:58 


164.12 


281.16 


9503 


2008-07-28 


17:54:40.2 


-29:36:60 


102.31 


275.21 


9504 


2008-08-02 


17:54:40.8 


-29:36:12 


125.42 


275.21 


9505 


2008-05-07 


17:54:37.5 


-29:34:47 


10.72 


82.22 


9854 


2008-07-27 


17:54:41.5 


-29:36:17 


22.78 


277.72 


9855 


2008-05-08 


17:54:37.5 


-29:34:47 


55.94 


82.22 


9892 


2008-07-31 


17:54:40.2 


-29:36:60 


65.79 


275.21 


9893 


2008-08-01 


17:54:40.8 


-29:36:12 


42.16 


275.21 


a Exposure time. 



52:00.0 40.0 20.0 17:51:00.0 50:40.0 

R.A. (J2000.0) 

Fig. 1. — Smoothed and exposure-corrected X-ray image (0.5- 
8.0 keV). All the ten data sets were st acked. The white circle 
at th e center shows the region studied by [Revnivtsev ct afj 1 120091 , 

120m 

sets using the Chandra Interacti ve Analysis of Obser- 
vations (CIAO) software package (jFruscione et al.ll2006l) 
version 4.2. We reprocessed the raw data to obtain an X- 
ray event list for each ObsID, in which we filtered events 
based on event grades, good time intervals, and energy 
bands (0.5-8.0 keV), and applied a background rejection 
algorithm to data taken with the very faint telemetry 
format 6 . We then merged the ten event files into one. 
Figure [T] shows the combined ACIS-I images corrected 
for the non-uniformity of the exposure across the study 
field. 

3. ANALYSIS 

3.1. Source Extraction 

We first extracted point source candidates using the 
wavdetect algorithm in the CIAO package. We set the 
significance threshold at 2.5 x 10~ 5 , implying that one 
false positive detection would be expected at every 4 x 
10 4 trials. As a result, 2,596 source candidates were 
found. The number of our source candidates is nearly 
the same with that by IRevnivtsev et al] (|2009D in the 
same region. 

For all the source candidates, we extracted source 
and background e vents using the ACIS Extract (AE; 
iBroos et al] I2010D package version 2009-12-01. The 
source events were extracted from a region containing 
^90% of photons based on the simulated point spread 
function of each source, while the background events were 
extracted locally from an annulus around each source. 

6 See http://cxc.harvard.edu/cal/Acis/Cal_prods/vfbkgrnd/index.html 
for detain 
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When two source extraction regions overlap with each 
other, each region was shrunk so that the overlap disap- 
peared. In the background regions, circles with a radius 
of 1.1 times the ^99% encircled energy radius of each 
source candidate were masked. All these processes are 
automated in the AE package, which is useful to extract 
sources in the studied region with a high source surface 
number density possibly with underlying extended emis- 
sion. 

To select significant point sources from the candidates, 
we examined their validity based on their photomet- 
ric significance (PS) and the probability of no source 
(Pb). The PS is defined as the background-subtracted 
source counts (C ne t) divided by its background counts 
normalized by the area. Pb is the probability that the 
source is attributable to a background fluctuation assum- 
ing the Poisson statistics. We recognized the source to 
be valid if they satisfy both two criteria: PS > 1.0 and 
P B < 1.0 x 10~ 2 . As a result, we obtained 2,002 valid 
point sources. Table [2] tabulates the photometric prop- 
erties of the first 20 sources, and the full source list is 
available in the on-line version. The X-ray sources can be 
referred following the International Astronomical Union 
(IAU) convention; e.g., CXOU J175044.88-292837.6 for 
the source sequence number 1 in Table [SJ 

3.2. Detection Completeness 

We estimated the detection completeness averaged 
over the entire area in the following method. We gen- 
erated events of 400 artificial sources using the marx 7 
event simulator at random positions across the image 
with a given flux. The events were merged with the ob- 
served events, from which point sources were detected 
using the same algorithm with the observed data. We de- 
rived the detection rate at each flux, which was changed 
from 10" 13 to 10 



~ 17 ergs cm 2 s 



in the total (0.5- 

8.0 keV) and the hard (2-8 keV) bands. We assumed 
that the artificial sources have the same spectrum with 
the best-fit model describing the stacked spectrum of all 
point sources (§ 14. 2j) . Figure [2] shows the detection rate 
as a function of flux. The rate is nearly complete (more 
than 97%) above ~10~ 15 ' 2 ergs cm -2 s^ 1 in the total 



band and ~10' 



-14.8 er g g cm 2 g 



in the hard band. 




10"'° 10"'° 10"" 10 ' 
Flux {ergs cm" 2 s 1 ) 



Fig. 2. — Detection rate as a function of flux in 0.5-8.0 keV 
(solid) and in 2—8 keV (dashed). 



3.3. Photometry 

Because most of the detected sources are faint, we used 
the photometric information to characterize their spec- 
tral hardness (colo r) and flux based on a quantile method 
(|Hong et al.ll2004l ). We derived the median energy (ME) 
for each source, which is a proxy for the spectral hard- 
ness. The distribution of the ME of all sources is shown 
in Figure [3] (a) . 




2 3 4 
Median Energy (keV) 



-16 -15 -14 

log F x (ergs cm" 2 s" 1 ) 



Fig. 3. — Distribution of (a) median energy and (b) photometric 
flux in the 0.5-8.0 keV band for all sources. 



We derived the flux of all sources using the photometric 
flux (Pph) defined as 



P ph (ergs cm" 2 s" 1 ) = 1.6 x 10 



_ 9 ME (keV) C net (s" 1 ) 
EA (cm 2 ) 



in which EA is the energy-averaged effective area for each 
source. We derived both the total (0.5-8.0 keV) and the 
hard (2-8 keV) band photometric flux. Some sources 
have no hard-band photometric flux because of their lack 
of photons in the band. The distribution of P p h is shown 
in Figure [3] (b) . In order to check the consistency with 
the flux derived by spectral fitting (§ 13. 5|) . we compared 
the flux estimates both by photometric and spectroscopic 
methods for 335 sources, for which spectral fitting result 
was available. The comparison shows that they are in a 
good agreement (Fig. Hereafter, we use the photo- 
metric flux for all sources. 




Flux by spectroscopy (ergs cm" 2 s" 1 } 

Fig. 4. — Relation between two X-ray flux estimates for 335 
sources, for which spectral fitting results are available. The verti- 
cal and horizontal axes represent the X-ray flux in the 0.5-8.0 keV 
estimated by the photometric and spectroscopic method, respec- 
tively. The dotted line represents the equal value between the two 
estimates. 



7 See |http: / /space.mit.edu/CXC/MARX/ 1 for details. 
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TABLE 2 

Chandra CATALOG: BASIC SOURCE PROPERTIES 





Source 




Position 




Extracted Counts Characteristics 




Seq 

# 
(1) 


CXOU J 
(2) 


R. A. 

(deg) 
(3) 


Decl. 

(deg) 
(4) 


Err 6 

(") (') 
(5) (6) 


C n et AC ae tC£ kg C net , hard PSFPS P B AnomVarEffExp ME Photo F x 

Frac (ks) (keV)(ergs cm -2 s —1 ) 
(7) (8) (9) (10) (11) (12) (13) (14) (15) (16) (17) (18) 


9i 92 
(19) (20) 



1175044.88- 


-292837.6267.687010 


-29 


.477113 0.411.6 


78.4 


23.5437.6 


29.2 0.89 


3.3 - 


3.9 




b 


486.5 


1.5 


1.6x10" 


-15 


-0.800.63 


2175045.10- 


-293431.4267.687940 


-29 


.575409 0.4 9.2 


38.4 


14.9163.6 


30.00.90 


2.5 - 


2.7 


g--- 




325.1 


3.0 


1.1x10" 


-15 


-0.811.26 


3175045.10- 


-293240.5267.687950 


-29 


.544603 0.3 9.7 


56.1 


19.1279.9 


17.70.89 


2.9 - 


3.3 


g--- 




535.5 


1.5 


2.4x10" 


-15 


-0.311.63 


4175045.51- 


-293222.4267.689630 


-29 


.539560 0.3 9.7 


46.4 


16.5201.6 


34.10.78 


2.7 - 


3.1 


g--- 




597.9 


2.9 


1.9x10" 


-15 


-0.331.57 


5175045.86- 


-293503.5267.691120 


-29 


.584322 0.3 9.1 


60.6 


16.5190.4 


52.70.90 


3.6 - 


5.0 


g--- 




417.5 


2.7 


3.0x10" 


- 15 


-0.391.67 


6175045.90- 


-293208.7267.691270 


-29 


.535760 0.3 9.7 


92.1 


20.3291.9 


24.20.85 


4.4< 


-5 


g--- 




639.9 


1.4 


1.5x10" 


-15 


-0.880.96 


7175046.03- 


-293539.1267.691800 


-29 


.594220 0.4 9.0 


40.2 


15.5179.8 


10.70.90 


2.5 - 


-2.7 


g--- 




387.2 


1.8 


1.3x10" 


-15 


-0.661.85 


8175046.07- 


-293405.8267.691960 


-29 


.568297 0.3 9.2 


69.0 


19.2270.0 


46.80.90 


3.5 - 


4.7 


g--- 




570.4 


2.4 


2.0x10" 


-15 


-0.461.40 


9175046.25- 


-292947.1267.692710 


-29 


.496437 0.310.7 


60.1 


17.6224.9 


51.20.69 


3.3 - 


4.3 




b 


678.2 


3.3 


2.5x10" 


-15 


-0.221.35 


10175046.62- 


-292935.0267.694250 


-29 


.493066 0.310.8 


58.6 


17.7229.4 


73.70.70 


3.2 - 


4.1 




b 


677.5 


5.2 


3.9x10" 


-15 


0.232.01 


11175047.08- 


-293334.4267.696170 


-29 


.559558 0.3 9.1 


48.4 


19.9316.6 


8.20.90 


2.4 - 


2.4 


g--- 




679.6 


1.2 


6.3x10" 


-16 


-0.961.12 


12175047.51- 


-293733.0267.697980 


-29 


.625834 0.3 9.0 


51.4 


17.4226.6 


42.10.90 


2.9 - 


3.3 


g--- 




519.2 


2.8 


2.0x10" 


-15 


-0.371.35 


13175048.28- 


-293250.7267.701180 


-29 


.547440 0.2 9.0148.7 


22.4321.3 


125.70.90 


6.5< 


-5 




a 


714.8 


3.5 


4.8x10" 


-15 


-0.181.34 


14175048.39- 


-292816.1267.701660 


-29 


.471139 0.411.3 


71.4 


20.2304.6 


47.4 0.90 


3.5 - 


4.5 


g--- 




386.5 


2.4 


3.1x10" 


-15 


-0.461.39 


15175048.54- 


-293937.5267.702280 


-29 


.660423 0.3 9.5102.3 


19.1235.7 


82.10.90 


5.2< 


-5 


g--- 




457.0 


3.0 


5.1x10" 


-15 


-0.291.45 


16175048.65- 


-293538.6267.702740 


-29 


.594083 0.2 8.5127.2 


18.6194.8 


79.10.80 


6.6< 


-5 


g--- 




734.4 


3.6 


4.8x10" 


-15 


-0.150.79 


17175048.86- 


-293736.9267.703620 


-29 


.626938 0.3 8.8 


73.8 


20.0296.2 


30.30.90 


3.6 - 


4.8 


g--- 




698.8 


1.9 


1.4x10" 


-15 


-0.650.84 


18175048.91- 


-292906.9267.703810 


-29 


.485275 0.310.6124.8 


25.5482.2 


45.00.90 


4.8< 


-5 




a 


668.0 


1.6 


1.9x10" 


-15 


-0.781.14 


19175048.95- 


-293547.9267.703980 


-29 


.5966610.2 8.5 


82.7 


16.4165.3 


51.40.76 


4.9< 


-5 




a 


752.1 


2.5 


2.1x10" 


-15 


-0.430.77 


20175049.07- 


-293803.6267.704500 


-29 


.634354 0.3 8.8 


82.8 


20.6312.2 


13.80.90 


3.9< 


-5 






702.2 


1.6 


1.3x10" 


-15 


-0.791.45 



NOTE. — Col.(l): X-ray catalog sequence number. £ 
standard deviation of the position error. Col. (6): Off- 
0.5-8 keV. Col. (10): Estimated net counts extracted in 
reduced PSF fraction (significantly below 90%) may indh 
the extracted counts (0.5-8.0 kcV) are solely from ba 
Col. (15): Variability characterization based o 
definitely variable {P K s < 0.005). No value is 
observed on-axis to obtain the reported nun 



orted by R. A. Col. (2): IAU designatior 
axis angle. Col. (7), (8): Estimated net 
the hard band (2.0-8.0 keV). Col. (11): 
:atc that the source is in a crowded rcgioi 



0.5-8.0 kcV band. Tablc[2] 
Col. (19), (20): qi and quantil 



blished : 
alu 



ckground. Col. (14): Source 
K-S statistic (0.5-8.0 kcV): a = 
eportcd for sources in chip gaps c 
icr of counts. Col. (17): Backgro 



. Col. (3), (4): R. A. and Decl. in the equinox .12000.0 
counts and errors extracted in 0.5-8.0 kcV. Col. (9): 
Fraction of the PSF (at 1.497 kcV) enclosed within the 
Col. (12): Photometric significance. Col. (13): Logi 



its 



entirety in 
icating the 



the 



-lcctn 



spectral shape 



jditic 
>f ca 



s: g = fractional time that source was on a detect 
vidence for variability (0.05 < P K s)\ b = possibly 
field edges. Col. (16): Effective exposure time: approximate tu 
:orrected median photon energy (0.5-8.0 kcV). Col. (18): Phot 
Astrophysical Journal. A portion is shown here for guidance rcge 
s (see text in 6 14.1.11 . 



!ol.(5)t Estimated 
kground counts in 
traction region. A 
ic probability that 
(FRACEXPO from mkarf) is < 0.9 
able (0.005 < P KS < 0.05); c = 



■ the 



flu 



uld ha 



the 



In order to examine time variability in the count rate, 
we applied the Kolmogorov-Smirnov (KS) test for the 
sources with net counts more than 100 (355 sources in 
total), which we call bright sources. The test examines 
the degree of non-uniformity in the distribution of photon 
arrival times against a uniform distribution (or constant 
flux) using the x 2 statistics. Sources were removed if 
they were in the chip gaps or at the field edge in at least 
one observation, because the telescope dithering causes 
artificial variability. 

The null hypothesis probability of the KS test was 
derived for each source in each of the ten observa- 
tions. We defined the merged probability (Pks) for each 
source as the smallest value among the ten probabilities. 
Based on this, we divided sources into three sets: (a) 
-Pks > 5 x 10~ 2 for non-variable sources (109 sources), 
(b) 5 x 10~ 3 < Pks < 5 x 10~ 2 for marginally variable 
sources (68 sources), and (c) Pks < 5 x 10~ 3 for variable 
sources (34 sources). 

By looking at the light curve of variable sources, we 
noticed that some sources show a typical flare-like flux 
variation of fast rise and slow decay. The KS test does 
not consider the profile of the variation in the light curve, 
so we picked up flare-like sources by setting a different 
criterion that the contrast of the maximum to minimum 
count rate is larger than 10 in at least one observation. 
In this manner, we selected nine flare-like sources, which 
are shown in Figure [5l 

3.5. Spectroscopy 



For the bright sources, we constructed the background- 
subtracted spectra and generated instrumental response 
files; i.e., redistribution matrix function (RMF) of the 
detector and auxiliary response file (ARF) of the tele- 
scope. We derived the best-fit parameters of the spectral 
models using the % 2 statistics. 

For the spectral models, we used an interstellar absorp- 
tion model (tbabs; Wilms et al. 2000) convolved with ei- 
ther of the two continuu m models: an optic ally-thin ther- 
mal plasma model (apec: ISmith et aT]l2001[ ) and a power- 
law model to represent the thermal and non-thermal 
spectra, respectively. Free parameters for the thermal 
model are the absorption column (-/Vh), the plasma tem- 
perature (k B T), and the flux (P x ) in the 0.5-8.0 keV 
band, while those for the power-law model are the ab- 
sorption column (-/Vh), the photon index (r), and the flux 
(Fx) in the 0.5-8.0 keV band. For t he plasma model, w e 
fixed the abundance to the value in Gu del et al.l (|2007l) . 
We regarded the fitting to be unsuccessful if the reduced 
\ 2 was larger than 1.5 or the best-fit values were unphysi- 
cal; i.e., T > 3 in the power-law fitting and k-&T > 15 keV 
in the apec fitting. For sources with unsuccessful fit- 
ting in both models, we also tried a plasma model with 
two different temperatures. For sources with successful 
fitting in both models, we adopted the result with the 
smaller reduced \ 2 value. 

As a result, 71 sources were successfully fitted with 
the one-temperature plasma model, 50 sources with the 
two-temperature plasma model, and 158 sources with 
the power-law model. Figure [B] shows the spectra and 
the best-fit models for all sources with more than 1000 
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Fig. 5. — Concatenated light curves for sources with flare-like variability. The light curves show variations in count rate (black histograms 
are binned and blue curves are adaptively smoothed). Time intervals be twee n observations are indicated by black dotted vertical lines. 
The top left labels show the sequence number in Table |21 and the group (§ 14.111 of each source. 



net counts. Tables [3] and [4] show the best- fit parameters 
for the one-temperature and power-law fitting of these 
sources, respectively. 

4. DISCUSSION 

We first attempt to group all the detected X-ray 
sources in § 14.11 based on the X-ray color-color diagram. 
After establishing these groups, we examine the X-ray 
properties of each group (§ I4.2p . We then discuss likely 
classes for sources in each group (§ I4.3|) and evaluate their 
contribution to the GRXE 

4.1. Grouping 

We cons t ructed an X-ray color-color diagram 
(Hon g et al.l [2004) using the quantiles (E25, E50, 
and £75) characterizing the spectral shape of each 
source. Here, E x (keV) is the energy below which x% 
of photons reside in the energy-sorted event list. E$q is 
equivalent to the median energy. We first normalized 
these parameters as 



E^ — Ey, 



E„ 



E„ 



(1) 



in which E m i n and i? ma x are 0.5 and 8 keV, respec- 
tively. We then took ratios of two parameters as q% = 
logio Q50/ (1 - Q50) and q 2 = 3Q25/Q75 for all the 
sources. Here, the q\ value indicates the degree of pho- 
ton spectrum being biased toward the higher (qi > 0) or 
lower (gi < 0) energy end (hard or soft spectra), and the 
92 value indicates the degree of photon spectrum being 
less (q2 > 1) or more (92 < 1) concentrated around the 
peak (broad or narrow spectra). 

Figure [7] shows the color-color diagram using q\ and q2 
of all the sources. In order to put the diagram into con- 
text, we simulated the quantiles of optically-thin thermal 
emission attenuated by interstellar photoelectric absorp- 
tion. We used the bremss model for the emission and 
the tbabs model for the absorption with varying tem- 
peratures and column densities. 

In the plot, sources are distributed from side to side 
with a main concentration around (qi, q 2 ) = (—0.8, 1.0). 
It appears that the scatter is extended in the upward 
and rightward directions from the main concentration. 
We examine if these branches have any physical basis by 
looking at count-stratified color-color diagrams (§ 14.1.1 1 
and the results of bright sources (§ 14.1.21 and § 14.1.3 ) 
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Fig. 6. — X-ray spectra of the sources wit h m ore than 1000 net counts in the order of the total counts. The top left labels show the 
sequence number in Tableland the group (§ [47TJ of each source. Thermal or non-thermal models are employed to fit the spectra. Grouped 
data are shown in the upper panels, which are over-plotted with the best-fit model convolved with the instrumental response. The lower 
panels show the residuals of the data to the fit. The best-fit parameters can be found in Tables \3\ and [4] 



before dividing all sources into groups by a statistical 
approach (§ I4.1.4[) . 

4.1.1. Count- Stratified Diagram 

The color-color diagram does not represent the flux 
of sources. Therefore, we constructed the diagram for 
different count ranges (Fig. [5J. The scatter plots ex- 
hibit quite a different morphology among the four count 
ranges. In the lowest count range (Fig. [Sh), the sources 
are scattered in the main concentration and the upward 
branch. As the count increases, the sources in the upward 
branch disappear and are gradually replaced with sources 
in the rightward branch (Fig. [5Jd) to form a V-shaped dis- 
tribution (Fig. [S]:). In the highest count range (Fig. 
the sources are only seen in the rightward branch. 

4.1.2. Variability of bright Sources 

In the variability analysis of the bright sources (jj 13.41) . 



some show flare-like variability. Figure [9] shows where 



they are located in the color-color diagram. We omitted 
the sources with less than 100 counts, as they are too 
weak to examine their variability; the faintest source with 
a statistically significant variability has a net count of 
100.5 (§ 13. 4p . It is interesting to find that all but one 
flare-like variable sources are found along the left arm of 
the V-shaped distribution, although the variable sources 
are equally found in both arms. 

4.1.3. Spectra of bright Sources 

In the spectral analysis of the bright sources (§ 13. 5| . 
some sources are described by thermal spectra and oth- 
ers by non-thermal spectra. For sources with more than 
1,000 counts (Fig. [6]), two sources are thermal and eight 
are non-thermal. Figure [TU] shows where they are located 
in the color-color diagram. The bimodality is clearly 
seen, in which thermal and non-thermal sources are ex- 
clusively found respectively in the left and right arm of 
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TABLE 3 

Thermal fittings for the X-ray sources with source counts over 1000 



Source 1 Spectral Fit 2 X-ray Flux 3 



Seq. CXOUJ Ct, M t Signif. log Nn m kT 2 F s F h F h , c F t F t>c x 2 /d.o.f. 

# (cm -2 ) (keV) (keV) (ergs cm -2 s _1 ) 

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) (13) 

287 175107.66-294037.3 3252.3 55.1 21.6±2;{£ 0.3±2;J 3.9+°;£ 3.6e-14 2.6e-14 2.6e-14 5.0e-14 5.2e-14 129.30/131 
1775 175148.91-293505.6 1779.1 40.7 21.6j][j;j£ 1.5^;J ... l.lc-14 8.1c-15 8.6c-15 1.9c-14 3.0c-14 103.52/ 78 

Note. — 1: For convenience, cols. (1)— (4) reproduce the source identification, net counts, and photometric significance data from Tabic [T1 
w hich are sorted by th e numb er of count s. 2: All fits used th e source model "tbabs^vapec" in XSPEC abundances frozen at the values relative 
to lAnders Grevess 1 J 1989ft . scaled to lWilms et alj J2000fl . using the tbabs absorption code in XSPEC. Cols. (5), (6), and (7) present the 
best-tit values for the extinction column density and plasma temperature parameters. Uncertainties represent 90% confidence intervals. More 
significant digits are used for uncertainties <0.1 in order to avoid large rounding errors; for consistency, the same number of significant digits 
is used for both lower and upper uncertainties. Uncertainties arc missing when XSPEC was unable to compute them or when their values 
were so large that the parameter is effectively unconstrained. Fits lacking uncertainties, fits with large uncertainties, and fits with frozen 
parameters should be viewed merely as splines to the data to obtain rough estimates of luminosities; the listed parameter values are not 
robust. 3: X-ray flux derived from the model spectrum arc presented in cols. (8)-(12): (s) soft band (0.5—2 keV); (h) hard band (2—8 keV); 
(t) total band (0.5—8 keV). Absorption-corrected fluxes are subscripted with a c. Cols. (8) and (12) arc omitted when log A^h > 22.5 cm~ 2 
since the soft band emission is essentially unmcasurable. 



TABLE 4 

Power-law fitting for the X-ray sources with source counts over 1000 



Source 1 



Spectral Fit 2 



X-ray Flux 3 



Seq. 
# 
(1) 



CXOU J 



(2) 



Ct.net Signif. log jV H T log N r 

(cm -2 ) 

(3) (4) (5) (6) (7) 



F s F h F h , c F t 

(ergs cm -2 s -1 ) 
(8) (9) (10) (11) 



Ft, c X 2 /d.o.f. 
(12) (13) 



1895 
438 
1149 
66 
33 
1831 
1252 
1869 



175154.71- 
175113.70- 
175131.68- 
175054.26- 
175051.16- 
175151.29- 
175134.06- 
175153.33- 



292806.1 
293110.1 
292957.0 
294327.4 
293419.5 
293310.3 
293103.9 
294245.0 



6440.6 
4325.4 
2908.3 
1600.2 
1469.3 
1447.3 
1390.3 
1202.7 



77.1 
64.4 
52.5 
46.2 
34.5 
36.3 
36.0 
30.0 



-, + (). 01 
-0.01 
-+0.04 
-0.04 
,+0.06 
-0.06 
7+0.04 
-0.05 
1+0.09 
-0.09 
■3+0.09 
-0.10 
,+0.08 
-0.08 
22 4+ 1 

zz -^-o.i 



22.0] 
22. 5j 
22.4j 
22.7j 
22.4j 
22.3j 
22. 1 H 



O + 0.1 

°-0.1 
r+0.1 
°-0.1 

Q + 0.1 

•°-0.1 
r+0.2 

'°-0.2 
7 +0.3 
'-0.2 

,DO_ 2 
fi +0.2 
D -0.2 

• U -0.3 



-4.26 
-4.53 
-4.87 
-4.87 
-4.92 
-5.51 
-5.18 
-5.30 



+0.06 
0.06 
+0.08 
0.08 
+0.10 
0.09 
+ 0.10 
0.09 
+0.2 
0.2 
+0.1 
0.1 
+0.10 
0.09 
+0.2 
0.2 



1.7e-14 1.4e-13 1.7e-13 1.6e-13 2.9e-13 300.71/244 

5.8e-15 l.lc-13 1.3e-13 1.2c-13 2.0e-13 196.80/203 

4.7e-15 7.8e-14 8.7e-14 8.2e-14 1.2e-13 137.31/142 

7.0e-16 8.0e-14 9.4e-14 8.0e-14 9.6e-14 128.96/110 

3.7e-15 3.7e-14 4.0e-14 4.1e-14 5.7e-14 88.60/ 71 

1.8e-15 4.8e-14 5.2e-14 5.0e-14 5.8e-14 113.47/ 78 

3.9e-15 2.5e-14 2.6e-14 2.9e-14 4.1e-14 75.74/ 74 

3.3e-15 5.0e-14 5.1e-14 5.3e-14 5.8e-14 60.59/ 58 



Note. — 1: For convenience, cols. (1)— (4) reproduce the source identification, net counts, and photometric significance data from Table [Tl 
which arc sorted by the number of counts. 2: All fits used the source model "tbabs*powcrlaw" in XSPEC. Cols. (5) and (6) present the best-fit 
values for the extinction column density and power law photon index parameters. Col. (7) presents the power law normalization for the model 
spectrum. Uncertainties represent 90% confidence intervals. More significant digits arc used for uncertainties <0.1 in order to avoid large rounding 
errors; for consistency, the same number of significant digits is used for both lower and upper uncertainties. Uncertainties arc missing when XSPEC 
was unable to compute them or when their values were so large that the parameter is effectively unconstrained. Fits lacking uncertainties, fits 
with large uncertainties, and fits with frozen parameters should be viewed merely as splines to the data to obtain rough estimates of luminosities; 
the listed parameter values arc unreliable. 3: X-ray Flux derived from the model spectrum, assuming a distance of 8.5 kpc, arc presented in cols. 
(8)— (12): (s) — soft band (0.5—2 keV); (h) hard band (2—8 keV); (t) total band (0.5—8 keV). Absorption-corrected fluxes are subscripted with a 
c. Cols. (8) and (12) are omitted when log ATh > 22.5 cm -2 since the soft band emission is essentially unmcasurable. 



the V-shaped distribution. 

4.1.4. Statistical Treatment 

From the discussion above, it is reasonable to as- 
sume that the sources consist of three representative 
groups that can be separated by the position in the di- 
agram: main concentration, upward branch, and right- 
ward branch. In order to define the boundaries in an 
objective manner, we employed the fc-means clustering 
algorithm using the mlpy script 8 . Given the number 
of groups, the algorithm determines a centroid of each 
group and the sources that belong to the group itera- 
tively, so that the sum of the distances from each source 
to the centroid becomes minimum. 

8 The details of t he script can be found in 
http:/ /mlpy. sourceforge.net/ 



Using this method, we divided the 2,002 sources into 
three groups (group A, B, and C) by their position in the 
diagram (Figure ITTj) . The resultant division was mostly 
consistent with its appearance; the group A for the right- 
ward branch, B for the main concentration, and C for 
the upward branch. Some exceptions are seen; e.g., two 
sources at the tip of the rightward branch were divided 
into the group C, rather than A, which is conceivable in 
such an algorithm. The two sources are quite faint com- 
prising only <0.5% of the total net count of the groups 
A and C, so they do not affect the composite spectrum 
of the groups discussed below. 

4.2. X-ray Properties of Each Group 
4.2.1. Fe K Feature 

We made composite spectra of the three groups A, B. 
and C using the CIAO tool combine_spectra. The tool 
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Fig. 7. — X-ray color-color diagram of all detected X-ray sources. 
The converted median energy is shown in the upper x-axis. We also 
display a grid of thermal plasma spectra attenuated by interstellar 
absorption with various parameter combinations. The labels indi- 
cate the values of the plasma temperature (k^T) for 0.5, 1,2,5, 10, 
and 30 keV and the extinction column (7V H ) for 10 20 5 , 10 21 , 10 21 5 , 
10 22 , and 10 22 ' 5 cm . In the simulation, background events were 
not included, so the grid should be considered only as a reference 
showing the relative parameter changes. 

makes background-subtracted spectra as well as count- 
weighted ARFs and RMFs of multiple sources. 

Figure [12] shows the enlarged view in the 4-8 keV 
band, in which the Fe K feature appears differently for 
each group. We characterized the feature by fitting the 
4-8 keV spectrum using a power-law plus a Gaussian 
model. The free parameters were the power-law index 
and normalization, and the Gaussian line center, width, 
and normalization. For the group C, the Gaussian com- 
ponent was not required, so we only constrained the up- 
per limit to the emission. The result is shown in Ta- 
ble [5] For the best-fit value of the Gaussian line cen- 
ter, we found no significant differences among different 
groups. 

The characterization of the Fe K feature was conducted 
also for the composite spectrum of all point sources and 
the total emission extracted from the entire field. The 
fractional contribution of the continuum (4-8 keV) and 
the Fe K band emission was derived separately for each 
group. The fraction of point sources among the en- 
tire emission is smaller than that reported (~80%) in 
iRevnivtsev et al.l fj2009f ) . which is presumably due to 
our lower sensitivity on average by using the entire field 
rather than the central part of the observations. 

The presence of the Fe K feature in the combined spec- 
trum of the group A sources (Fig. [T2k ) seemingly contra- 
dicts the result of the spectral fitting of individual bright 
sources (§ I3.5|) . because most of the bright sources in 
the group A hardly show signatures of the Fe K feature. 
Among the brightest ten sources shown in Figure El eight 
are group A sources (source numbers 1895, 438, 1149, 66, 
33, 1831, 1252, and 1869), all of which are fitted with a 
power-law model (Table 0]) . This strongly suggests that 
the presence of the Fe K feature depends on the source 
flux within the group A. No such apparent contradiction 
was found in the other groups. 

In order to investigate the flux dependence for the Fe K 
contribution, we constructed spectra stacking flux-sorted 
sources cumulatively from the brightest source toward 
the fainter end by increasing the source number by 1 (1- 
20), by 10 (20-40), and by 20 (40-) at each step. In 
other words, the composite spectrum was made for the 



brightest 1, 2, 20, 30, 40, 60, 80, ... sources. We 
then applied a power-law plus a Gaussian model for the 
spectrum generated at each step to derive the EW of the 
Fe K feature. Figure [T2] (a) shows the EW value against 
the decreasing flux. At the brightest end (>2xl0~ 14 
ergs cm~ 2 s _1 ), the EW is consistent with being zero, 
which gradually increases as the threshold flux decreases 
and eventually levels off. This suggests that different 
classes of sources are dominant in the brighter and fainter 
end of the flux and their ratio gradually changes along 
the flux. 

4.2.2. Global Spectral Fitting 

We performed global spectral fitting using the entire 
energy band (0.5-8.0 keV) for the composite spectrum 
for each group (Fig. 1111) in a similar manner for the 
individual bright sources (§ 13.41) . For the thin-thermal 
pl asma model, we fix ed the abundance to that adopted 
in lGiidel et atl (fl999l) . 

Group A — For the group A spectrum, neither a power- 
law nor a thin-thermal plasma model did not reproduce 
the spectrum well respectively because of the excess emis- 
sion at 6.7 keV or a flatness of the continuum. In fact, 
these two are a signature respectively of a thin-thermal 
and a power-law model, so we fitted the spectrum with a 
combination of the two, which was successful. The best- 
fit model are given in Figure [TJJ whereas the best-fit 
parameters are in Table |H1 

The EW of the Fe K feature becomes larger as the flux 
decreases (§ l4.2.ip . which suggests that the thermal com- 
ponent becomes stronger against the non-thermal com- 
ponent as the flux decreases. In order to investigate 
this trend, we fitted the cumulatively stacked spectra 
(§ 14.2. ip with a power-law plus an optically-thin thermal 
plasma model. The free parameters are the normaliza- 
tion of the two components, and the other parameters 
were fixed to the values obtained in the fitting of all group 
A source (Table [5]) . As is expected, the flux ratio of the 
two components -F p0 wer-iaw/-Fthormai starts to decrease 
at ^-TO -14 ergs cm~ 2 s _1 fFig.lTB'b'). which is consistent 
with the EW trend fFig.lTBlO. 

Group B & C — For the group B spectrum, several emis- 
sion lines are seen, including the 6.7 keV emission from 
Fe XXV and 2.5 keV from S XV. This set of emission 
lines indicates a multiple-temperature plasma, and in- 
deed the spectrum was reproduced well with two thin- 
thermal plasma components, but not with one compo- 
nent. The free parameters are the plasma temperature, 
metallicity, and absorption column. The metallicity was 
the same among all the metals and was constrained to 
be the same between the two plasma components. 

For the group C spectrum, several emission lines are 
also seen. Unlike the group B sources, the Fe XXV emis- 
sion at 6.7 keV is absent. We fitted the spectrum using 
the same model with group B and obtained the best-fit 
parameters in Table [6] 

4.2.3. The log ./V-log S Curve 

We made the log TV-log S curve separately for all point 
sources and those for each group (solid histograms in 
Fig. 1 15[) . The curves were constructed in the hard band. 
For the group A curve, we divided the cumulative number 
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Fig. 8. — X-ray color-color diagram for different count ranges. The symbols are the same with Figure I7l 



TABLE 5 

Best-fit parameters of the power-law + Gaussian model for the Fe K FEATURE a 



Group 


JT b 

-C'gaii 

(keV) 


EW g au C 

(eV) 


Fpow/107 13d 

(ergs cm -2 s -1 ) 


f (PS) e 
Jpow 


.(Total) f 
Jpow 


F gau /10- 13 s 
(ergs cm -2 s —1 ) 


r(PS) h 
J gau 


..(Total) i 
/gau 


A 
B 

C 


6 69+ 06 
°- DM -0.05 

6 70+ - 08 
°- ,u -0.08 


2681^ 


12 02 +a58 

b-oi1S:S 

o.68lg;g 


U - o8 -0.03 
n 9s +0.04 
U - zS -0.03 
04+ 01 


u -^-0.02 
n 1s +0.02 
U - 18 -0.02 
02+ 01 


<0.04 


62+ 08 

U ' DZ -0.08 

0.3210-22 


26+0- 03 
u - ZD -0.03 
iq+0-09 


All PS 

Total emission 


6 69 +aub 
D - Da -0.04 
6 74+0.05 


oni+llb 

oz -122 

58iii°; 


97 99 +0.54 
zi.zz_ 54 






i ro+U.bU 
1 - t,8 -0.51 
n on+0.50 
d - 8z -0.51 







a The best-fit value and a ler statistical uncertainty are given. For those with no best-fit values, only the 1 a upper limit is given, 
k Center energy of the Gaussian component. 
c Equivalent width of the Gaussian component. 

Flux of the power-law component in the 4—8 keV band. 
c Fractional flux of the power-law component among all point sources. 
^ Fractional flux of the power-law component among total GRXE spectrum. 
s Flux of the Gaussian component. 

_ Fractional flux of the Gaussian component among all point sources. 
] Fraction flux of the Gaussian component among total GRXE spectrum. 



TABLE 6 

Best-fit parameters for global spectral model in 0.5-8.0 keV 3 



Group 


JV H (1)b 
(10 22 cm" 2 ) 


(keV) 


AT H (2)d 
(10 22 cm" 2 ) 


fc B T( 2 ) c 
(keV) 


z f rs 


N W /Ar (2) h 

1 v apcc / 1 Y apcc 


X 2 /d.o.f. i 


A 
B 

C 


1 09+0- 39 

o.75±8:gg 

0.701°;}? 


6 65+ 3 ' 24 
0-00-3. 03 

74+0-54 
U -' 4 -0.45 
n 7 o + 0.04 
U -'°-0.03 


2-46lg;I 

0.80l°; 22 

04+ 05 
u - u4 -0.04 


7 07+I.86 
'•o'_4.84 

4 50+0- 65 
4 - ou -0.35 


O7+0-36 1 2q +0.18 

u - M '-0.32 i - zy -0.40 

qq+0- 33 
u.yy_ 29 

15+ 018 


47+0- 18 

6 70+0- 13 
O- 'O-0.20 


205.36/504 
98.68/103 
85.86/102 



a The best-fit value and a lcr statistical uncertainty are given. 

k Interstellar extinction column density for the first component (the lower temperature component for the two- temperature model). 
c Plasma temperature for the first component (the lower temperature component for the two-temperature model). 

Interstellar extinction column density for the second component (the higher temperature component for the two-temperature model). 
e Plasma temperature for the second component (the higher temperature component for the two-temperature model). 
^ Metal abundance relative to the solar value for the thermal component. 
s Photon index for the power-law model as the second component. 

Normalization ratio between the first (the lower temperature) and and the second (the higher temperature) component for the 
two-temperature model. 



1 The goodness of the fit. 

into two, representing the power-law and plasma sources 
by using the smoothed flux ratio curve as a function of 
flux (Fig. 113b). We did not correct the log TV-log S curve 
because we do not discuss the absolute value of the ver- 
tical axis of the log TV-log 5 curve. 

It is expected that a certain fraction of point sources 
are extra-galactic sources, which we represent with the 
logiV-l og5 curve construc ted in a Chandra deep field 
(CDF; IRosati et al.1 l2QQ2h . Since the extra-galactic 



sources in our study field are more absorbed than those 
in the CDF at a high Galactic latitude, we took into ac- 
count the flux reduction due to the Galactic absorption 
with a column density of 10 22 cm -2 assu ming a typi 



cal s pectral shape of sources in the CDF (Tozzi ct al. 
2006). The additional Galactic column density is a typi- 
cal best-fit value among the sources with spectral fitting 
(Table El). 
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Median Energy (keV) 

2 3 4 5 




-1.4-1.2-1.0-0.8-0.6-0.4-0.2 0.0 0.2 0.4 
q^log^Q^/fl-Q^) 



Fig. 9. — X-ray color-color diagram of the bright sources (Cnet 
> 100 counts). Variable sources are shown with magenta circles, 
among which flare-like variable sources are filled. 



Median Energy (keV) 

2 3 4 5 




-1.4-1.2-1.0-0.8-0.6-0.4-0.2 0.0 0.2 0.4 
q^logroOso'd-Qso) 



Fig. 10. — X-ray color-color diagram of the bright sources (Cnet 
> 100 counts). Sources with a thermal spectrum (Fig. [6)l are 
shown with red circles, while those with a non-thermal spectrum 
are shown with blue circles. 



Median Energy (keV) 

2 3 4 5 




-1.4-1.2-1.0-0.8 -0.6-0.4-0.2 0.0 0.2 0.4 
qi=l09ioQ50 / ( 1 -°5o) 



Fig. 11. — X-ray color-color diagrams of all sources color-coded 
depending on their statistically-defined classes. The other symbols 
follow Figure [71 




Fig. 12. — Fe K feature in the combined spectra of point sources 
in the group (a) A, (b) B, and (c) C. The spectra are normalized 
at 4 keV. 




F x (ergs cm" 2 s" 1 ) 



Fig. 13.— (a) EW of the Fe K feature and (b) the 2-8 keV 
flux ratio of the power-law (-F p0 wer— law) versus thermal (-^thermal) 
components against the decreasing flux in 2-8 keV above which the 
cumulative combined spectra were constructed. The lcr statistical 
uncertainty is shown for each data point. 



4.3. Likely Populations 

We now discuss likely classes of the X-ray point sources 
in each group based on the results presented above. The 
transition from one class to another is continuous along 
the color and flux, so the groups are naturally a mixture 
of sources of different classes. We thus discuss the likely 
classes of sources that constitute the majority of each 
group. 

First, we consider that the group A sources are mostly 
a mixture of AGNs and WD binaries, each responsi- 
ble for the power-law and thermal plasma component 
in the composite spectrum (Fig. nib). For the power- 
law component of the group A, the \ogN-\ogS curve 
(Fig. [T5)) matches quite well to the compl eteness limit 
(Fig HI) with that obtained in the CDF (jRosati et al.l 
2002). In the spectral fitting, the power- law component 
has a spectral index of 1.29^40 (Table [S]), which is simi- 
lar to the typical spectral form of AGNs with a power- law 
index of about 1.63±0.13 in the X-r ay band regardless 
of their types and lumi nosities (e.g.. iRosati et all 120021 : 
ICharles fc Seward [1991 . 

For the thermal component of the composite spec- 
trum (Fig. H3k) of the group A, a strong Fe K feature 
and a 7 keV temperature plasma (Table [6]) are seen. 
Both of these feat ures are observationa l characteristics 
of magnetic CVs (|Ezuka fc Ishidalll999h . Other classes 
of white dwarf (WD) binaries, such as dwarf novae, pre- 
CVs, and symbiotic stars, may also be considered for 
the likely classes. Pre-CVs are poorly recognized class 
of sources, which are detached binaries of a WD and a 
late-type star, unlike conventional CVs that are semi- 
detached systems. Som e of them show strong Fe K emis- 
sion in the hard X-rays (Matranga et al. 2012). Likewise, 
some symbiotic star s are also know n to show strong Fe 
K emi ssion feature (|Eze et al.ll2012f ). In fact. iMorihanal 
(2012) showed near-infrared spectra of selected X-ray 
sources presented in this paper, in which some thermal 
A sources do not exhibit the Br7 emission that i s typ- 
ical for the conventional CVs (jDhillon et al.l 119971 ). We 
thus consider that pre- CVs and sym biotic stars, most of 
which do not show the Biry em ission (|Howell et al.ll2010t 
lSchmidt~fc Mikolaiewska 2003|), also account for at least 
some fraction of the thermal source population in the 




Fig. 14. — Composite spectra and the best-fit global model of the group (a) A, (b) B, and (c) C. The symbols follow Figure [6] 




Fig. 15. — Cumulative source surface number density N (deg 2 ) 
against the flux S ergs cm" 2 s" 1 in the hard band (2—8 keV) sep- 
arately for (a) all point sources (green) and those belonging to the 
group A (blue), B (black), and C (magenta), and (b) those decom- 
posed into power-law (A powcr _i aw ; red) and thermal (At^ermai; 
cyan) sources. The curve for extra-galactic sources with an ab- 
sorption of ~10 22 cm" 2 is shown w ith the gray dotted curve using 
the CDF data URosati et al.ll2002l) in (a). The slope of -1.5 for 
spatially uniform distribution is shown at the top right. All curves 
are not corrected for the detection completeness. The nearly (97%) 
completeness limit is shown with the vertical dashed line. 

group A. 

The group B and C sources are Galactic sources with 
a soft thermal spectrum, and we consider that most of 
them are likely to be X-ray active stars. The composite 
spectra of the two groups were fitted with two plasma 
components (Fig. [Tj] and Table . The lower temper- 
ature component and its absorption column are quite 
similar between the two, indicating that the difference 
between the B and C are not their typical distances. On 
the contrary, the higher temperature component is sig- 
nificantly higher for the group B; it is is high enough to 
produce a strong Fe K feature at 6.7 keV, which is ab- 
sent in the C spectrum. We speculate that the difference 
between B and C is that most C sources represent ac- 
tive binary stars in the quiescence, while most B sources 
represent those during flares. Although most flare-like 
variable sources are indeed group B sources (Fig. this 
can be a statistical bias because the group C sources are 
systematically fainter than group B sources (Fig. [8]). 

4.4. Comparison with Previous Studies and 
Contribution to GRXE 

We have constructed a log iV-log S 1 curve in the 2- 
8 keV band down to the flux ~1 -16 ergs cm -2 s _1 using 
the CBF data (§ 14. 2\ . Several major groups of sources 
account for different fractions of the GXRE continuum 
and Fe Ka emission in different flux ranges. 

For the continuum emission, low-mass X-ray binaries 
dominate in the flux range above «10" 12 ergs cm" 2 s" 1 
in the log iV-log S curve, as was found in previous studies 



(e.g.. lHands et al.l l2004). These bright Galactic sources 
saturates, and background AGNs emerge as the most 
dominant class of sources in the next flux range of 
?sl0" 12 -10" 14 ergs cm" 2 s" 1 . The dominance of AGNs 
over Galactic population in this flux r a nge was also 
discus se d bv lEbisawa et al l (120011 120051): lHands et al.1 
(|2004[ k iRevnivtsev et all (|2009f ): iMotch et al.1 (|2010T ). 
Below ?sl0" 14 ergs cm" 2 s" 1 , Galactic sources start to 
dominate the \ogN-\ogS curve again. We decomposed 
the log TV-log S curve of the group A sources into non- 
thermal (AGNs in our speculation) and plasma sources 
(WD binaries) by utilizing the increasing EW of the Fe 
Ka emission in the composite spectra along the decreas- 
ing flux. The group A thermal sources (WD binaries) 
and the group B sources (X-ray active stars) are the ma- 
jor contributers in this flux range. These two Galac- 
tic class es are also proposed to be dominant i n this flux 
range bv IRevnivtsev et ail (|2009( ) : IHongl (|2012l ). The for- 
mer study claimed that the X-ray active stars account 
for more than a half of the integrated flux, while the 
latter study claimed the opposite. In our study, the inte- 
grated number of sources of the group A thermal sources 
(WD binaries) are taken over by that of the group B 
sources (X-ray active stars) at ~3xl0 -15 ergs cm" 2 s" 1 
(Fig. El. 

For the Fe Ka emission, we found that the group A 
thermal sources (WD binaries) and the group B sources 
(X-ray active stars) are the two major contributers (Ta- 
ble H), despite their much smaller contributions to the 
continuum emission in comparison the group A non- 
thermal sources (AGNs). If we assume that all the Fe 
Ka emission of the group A sources (Fig. IT4"k ) is from 
its thermal component, the fractional contribution to the 
Fe Ka emission by the group A thermal (WD binaries) 
and group B sources (X-ray active stars) is about 2:1 
(Table EJ. 

5. SUMMARY 

We presented the result of coherent photometric and 
spectroscopic analysis of 2,002 X-ray point sources us- 
ing the entire data set of the deep Chandra bulge field 
observations, based on which we discussed X-ray source 
population in the hard X-ray band in the flux range be- 
tween wlO -13 ergs cm" 2 s and «10 _16 ergs cm" 2 s" 1 
contributing to the GRXE. The main results are as fol- 
lows: 



1. Based on the X-ray photometric information, we 
phenomenologically classified all the sources into 
three groups A, B, and C. The group A sources are 
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the mixture of two different classes of a thermal and 
non-thermal spectrum. From their X-ray proper- 
ties, we speculated that the most dominant class 
of these groups are (1) background AGNs (group 
A non-thermal), (2) WD binaries (group A ther- 
mal), (3) X-ray active stars during the flare (group 

B) and (4) X-ray active stars in quiescence (group 

C) . 

2. We made a composite X-ray spectrum of each pop- 
ulation and derived their fractional contributions 
to the GRXE separately for the hard-band contin- 
uum and the Fe K emission line feature. We found 
that the group A thermal and group B sources ac- 
count for most of the Fe K emission. 

3. We constructed the \ogN-logS curves 
for each group in the flux range between 
«10~ 13 ergs cm~ 2 s _1 and «1CP 16 ergs cm~ 2 s _1 
in the 2-8 keV band. We found that extra- 
Galactic sources dominate the flux range down to 
ss 5 x 10~ 15 ergs cm~ 2 s _1 , which is taken over by 
the Galactic populations below this flux. For the 
Galactic sources, WD binaries such as magnetic 
and non-magnetic CVs, pre-CVs, and symbiotic 
stars, account for a larger fraction than X-ray 
active stars in the studied flux ranges. 
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